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Hie work described in this report wee perforaed for the Auxiliary 
Propuleion Branch of Che Propulsion Division, Structures and Propulsion < 

Laboratory, NASA George C* Marshall Space Plight Center^- under*^Contract | 

NAS8*33974. The NASA technical monitors were Lee W. Jones and T.D. McCay. 1 
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This report represents Part II of < the Pinal Report. Part 1 deals with 
the evaluation of alternate laser devices as well as alternate propellants 
and energy coupling schemes. 
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ABSTRACT 

Ttie present etetue of theoretical aodela for the laser heated thruster 
is reviewed. It is concluded that exiating aodela neither agree with each 
' Owher nor with the limited experimental data available. The requirements 

for an improved laser heated thruster theoretical model era discussed. The 
application of a time-dependent finite-difference Naviar-Stokea equation 
solution to the laser heated thruster problem is described, along with a 
simple closed form solution which was developed in order to gain insight 
into tiiw difficulties encountered in the pursuit Of the numerical solution. 
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1. LASER HEATED PROPULSION MECHANISM 


The concept of laser heated propulsion is based on using a laser beam 
to remotely heat a working medium in a rocket thruster. Such a scheme 
appears to be attractive because it might be possible to generate propellant 
gas temperatures which are much higher than those possible in chemical pro- 
pulsion devices. In particular, in conjunction with propellants of low 
molecular weight, a very high thrust chamber temperature will result in a 
specific impulse much higher than that achievable in combustion driven 
rocket motors. Even through the concept ir simple, the sttendsnt high 
Cenperacure presents difficult engineering problems. 

The use of continuous wave lasers for laser heated propulsion can be 
divided into s very high specific impulse, high temperature regime and a 
high specific impulse, awdium temperature regime. While the first regime 
would be realised by utilizing hydrogen as the propellant and a laser 
supported absorption wave generated and sustained by inverse bremsstrahlung 
absorption, in the medium temperature regime laser energy would be coupled 
into a suitable propellant by molecular or aerosol absorption. Details of 
these and other schesMS are discussed in Part I of this report. 

The work described here is mainly concerned with an assessment of laser 
propulsion fluid mechanics modeling and simulation requirements for the very 
high specific impulse, high temperature regime. Here the essential feature 
IS that the propellant gas, which is normally a transparent, non-conducting 
medium, becomes ionised by the laser beam and then strongly absorbs energy 
from the beam (Ref. 1). 

If a few free electrons are present in the focal volume, they can gain 
sufficient energy from the beam to produce further ionisation via collisions 


with neutral gas particles. Breakdown will occur if the electron density 
reaches a critical value despite losses due to diffusion., ittachment, 
recombination and other nechanisas. 

Oiioi- <1 plnsDui is.ts boon fumed, it eaiCh ruiliation whoso wavelength dis- 
tribution depends on the temperature and the density. Mo»t of Che energy 
emitted is in the far ultraviolet and is believed to be more readily 
absoibed in Che surrounding cold gas racher chan in Che plasma icself. 

Thus, a layer of gas ouCside Che plasma, although transparent to laser 
radiation, is heated by the plasma radiation, and upon reaching a suifi- 
ciently high temperature, will be ionised to such ar. extent that it will 
also become absorbing for the laser beam. This layor will then be further 
heated rapidly until its temperature beco’ses so high that it in turn will 
become transparent to ultraviolet plasma radiation. By this time a new 
layer of plasaui nearer the laser will have become absorbing, so the boundary 
of the plasms, called an absorption wave, will m>.>ve toward the laeer. If 
this occurs in a flow of propellant gas, the tlow velocity of which is equal 
in iiKigntlude hnt opposite in dii.-vtion to the motion ol the obeorptiun wove, 
we have a stationary propellant plasma which can function as the heating 
mechanism for a propulsion device. 

Various laser heated propulsion devices have been conceived and 
proposed. Since it would be difficult to discuss the modeling and simula- 
tion requireownts in complete generality, one particular configuration for a 
laser heated thruster is shown in Fig. I (schematically). This is the 
arrangement used in the aubsequent discussion of sKideling and simulation 
requirements. In particular, we are interested in describing or modeling 
the interaction of the initic’ly cold propellant gas (H^) as it flows 
through and around the laser supported plasma region, then mixes and exits 
through the subsonic-supersonic nottle, producing thrust in the process. 
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Focused Lskser Beam 



F ijj 1 - Laser Heated Propulsion Device (Schematic) 
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Kor llitf .icttial of such a device il is most inportanl Co predict 

the propellant velocity required to achieve a stable plaaaa as well aa Che 
cooling requirenents for the plasma and mixing chamber walls. Hydrogen is 
of particular interest as a propellant because it promises Che highest 
specific impulse. 


7. REVIEW OF EXISTING MODELS 


2.1 ABSORPTION WAVE MODEL 

Although nuch has been achieved in modeling the physics of Che laser 
heated thruster, much must be done toward obtaining agreement with experi- 
ment if these mathemaCioil models are Co serve as useful Cools in the design 
of laser heated propulsion devices. Most of Che existing models deal with 
absorption waves in air and have numerous simplifying assumptions in coimon. 
Verv little work so far has been done with hydrogen as the propellant. 
Experimental data for absorption waves in hydrogen are also very sparse. 

The first model Co appear in the literature was Chat of Raiser (Ref. 

2). His model is based on the assumptions of one-dimensional flow at 
constant pressure with heat conduction considered to be the principal heat 
transfer mechanism. The plasaia is considered to be optically thin to its 
self-radiance, and in thermal equilibrium. This allows ths electron density 
Co be determined from the Saha equation. Also, Che absorption coefficient 
for the absorption of laser radiation via inverse bremsstrahlung as wall as 
an approximate function for the overall radiation losses can be specified as 
functions of temperature and pressure. 

Finally, in order to arrive at an analytical solution, the absorption 
coefficient as well as the ratio of specific heat to thermal conduction are 
asitumed constant. The neglectiun of radiation losses limits the applic- 
ability of the solution to small plasma volumes (diameters less than a few 
millimeters at p ■ I Stm). 

With the goal Co improve upon Che analysis by Raiser, Jackson and 
Nielsen (Ref. 3} presented a model which includes the mechanism of radiative 
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transfer into Che ambient gas (air). Apart from this, other assumptions 
made are basically the same as those invoked by Raiser. Jackson and Nielsen 
derive an analytic solution, the evaluation of which is performed by iterat- 
ing cn Che value of the radiative loss function. The authors point out 
seveial shortcomings of their analytical solution: (1) interferograms taken 

of absorption waves show clear evidence of radial flow, rendering Che prob- 
lem two-dimensional, and (2) their calculated propagation velocity of Che 
absorption wave as a function of incident laser intensity does not agree 
with experiment. 

Shortly following the analysis of Jackson and Nielsen, batceh and 
Keefer (Ref. 4) presented a two-dimensional generalization of Raiser 's 
analysis for the subsonic propagation of laser spaiks in air. To be pre- 
cise, BatCeh and Keefer combined a one-dimensional flow velocity description 
with a Cwo-diMnsional temperature field description, retaining Raiser's 
assumption of ordinary thermal conduction as the principal propagation 
mechanism, with plasoui radiation serving only as an energy loss mechanism. 
The importance of this two- limens ions 1 temperature field model is chat it 
provides insight into the effect of boundaries on spark propagation in a 
channel. Other assumptions such as constant pressure, constant ratio of 
specific heat to thermal conductivity are also retained. The absorption 
coefficient was assumed to be piecewise constant, zero in front of the 
absorption wave, and non-s«ro within the plasma region. An analytic solu- 
tion in terms of Bessel functions was obtained to evaluate the relationship 
between laser intensity and spark propagation velocity. 

BatCeh and Keefer found that their calculated prop>ii{ation speeds are 
roughly an order of niag(<icude lower than those obtained experimentally, 
while the maximum temperature in the wave as calculated matcheo that of the 
experiment. The authors suggest that this is due to considering thermal 
conduction only as the primary propagation mechanism. They essentially 
argue th-it plasma radiation, causing ionization and additional laser energy 
absorption via inverse bremsatrahlung in front of the absorption wave must 
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be included in Che model Co achieve realiscic predicCiona of Che propagacion 
velocicy as a funcCion of laser intensicy. 

The i.exc model, thac of Kenp and Root (Ref. 3), in fact does include 
Che effecC of forward plasma radiacion in a one-dimenaional numerical 
model. This model assumes one-dioiensional flow ac conscanc pressure in a 
conscant area channel. Since che equacions are solved by numerical incegra- 
tion. It IS possible to treac the thermodynamic properties as well as the 
various transport phenomena (absorption, thermal conduction and radiation) 
in a much more detailed manner chan in the previously discussed analytical 
solutions. The solution of Kemp and Root's formulation requires an itera- 
tion on the mass flux for specified laser intensity. Only one particular 
mass flux value will yield a physically reasonable temperature distribution 
in Che absorption wave (saddle-point singularity). Calculations are pre- 
sented for pressures of 1, 3, 10 and .10 atm, and for laser powers of 10 kW 
and 3 HW. The authors observe a very steep rise in temperature at the 
leading edge of the absorption wave and cautio.i that for such conditions the 
"radiacion conduction approximation," used to consider the heating effect of 
plasma radiation on Che cold gas in front of che wave, is not really valid. 
However, since energy is properly conserved in rhe model, they suggest that 
the temperature profile in the wave proper should be qualitatively correct. 
No experimental data for absorption waves in hydrogen were available at the 
time to confirm the theoretical predictions. Also, the authors realise that 
che one-dimensional nature of the model presents serious limitations on the 
validity of the results. It is pointed out that thermal conduction in the 
radial direction becomes important for small plasma regiona ar low pressure. 
A very obvious two-dimensional phenomenon is the focused laser beam. In 
this case, the intensity will not only change due to absorption but also due 
to varying cross-sectional area. A convergent beam behaves as a higher 
intensity beam - relative to a collimated beam of the same initial intensity 
- and requires a higher propellant flow velocity to keep the "plasmatron" 
stationary. Finally, there is the transverse velocity component. It is 
suggested that its effects should be less in a confined flow - such as in a 
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thruster - than in the unconfined case. Since neither theoretical nor ex- 
perimental results are available, it is difficult to assess the importance 
of two-dimensional flow for absorption waves in a channel. 

In concluding this brief discussion, it appears that the modeling of 
absorption waves is still in a stage of infancy. In Fig. 2 (adopted from 
Ref. 3) we have sumnarized results for sbsorption waves in air from the work 
reviewed here. Some experimental results are also shown. No attempt has 
been made to systemat ical ly correlate these results. It is evident, How- 
ever, that theoretical predictions are in considerable disagreement with the 
experimental data (Refs. b,7) and with each other, even for similar assump- 
tions in the theoretical models. 

2.2 THRUSTER AND NOZZLE FLOW MODELS 

In order to produce thrust, the heated gas exiting from the absorption 
wave must be accelerated through a supersonic noztle which basically 
converts thermal energy into kinetic energy. Just as for the absorption 
wave, the modeling of this flow can bu approached at different levels of 
sophistication. 

The first of such models to be discussed here is that of Kemp and Root 
(Ref. 5). It is a one-dimensional model, obtained by explicitly including 
the area term and retaining the momentum equation in the set of one- 
dimensional equations used for the one-dirnsnsional absorption wave model by 
the same authors. The results of the absorption wave calculations serve as 
initial values for the thruster problem, i.e., the soluCi ns are joined at 
an appropriate location downstream of the wave. The model described in Ret. 
5 represents an "inverse method." In Ref. 5, rather than specifying the 
thruster and nozzle contour, the velocity is specified as a function ot flow 
distance. While this or variations of this procedure are coinnon in order to 
avoid the sonic singularity at the tnrost, it means that many iteractons mz) 
have to be performed until a reasonable or specified nozzle contour is 
obtained . 
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While Che one-dimeneional chruater flow model haa the advantage of 
simplicity, it cannot properly account for lateral heat condution and radia- 
tion Co Che walla of the thruaC chamber. As pointed out in Ref. 3, a pre- 
cise evaluation of the transverse heat loss due to thermal conduction and 
radiation requires a two-dimensional model. Such a model was presented by 
Che same author (et al.) in Ref. 8. This model is based on Che boundary 
layer equations and thus includes radial gradients. Added terms for laser 
energy absorption, eqjilibrium chemistry and plasma radiation losses 
complete the equations. Although this model is an improvement over the 
previously mentioned o.ie-d imensional model, it does not contain any terms 
for forward energy transport by rauiation or conduction, and therefore 
cannot be applied to model the absorption wave. A method to join the one- 
dimenaional absorption wave calculations to the two-dimensional thruster 
calculations is described in Ref. 9. According to the authors, the two- 
dimensional thruster code needs a number of refinements to ouke it generally 
useful. Required are a general tsation of the radiation model, more exper- 
ii'ncc or a rationale to choose the axial pressure gradient to obtain a 
desired flow channel contour, and a variable coordinate grid structure to 
concentrate grid points in areas of steep gradients. 

2.3 THERMODYNAMIC AND TRANSPORT PROPERTIES 

As compared with other fluid flow problems, the description and 
analysis of the flow in a laser-heated thruster is severely complicated by 
the large temperature rang..s over which thermodynamic and transport phe- 
nomena have to be considered. As a consequence, the early investigations, 
in particular those that produced analytical solutiona, had to employ many 
simplifying assumptions (Kets. 2, 3 and <*) , Ultile the qualitative nature of 
til.' ri’siilts from these .inalyses is useful in furthering an understanding ot 
th<> phenomena involved, the pursuit of quantitative results probably denumds 
mimertcal solutions and imicli more detailed thermodynamic and transport prop- 
erties, such as used by Kemp et al. (Refs. 3, 8 and 9). In fact, s substan- 
tial portion of the work by Kemp et al. deals with the development, evalu- 
ation and assembly of supporting data .such as laser absorption coefficients. 
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as well as thermodynamic and transport properties (conduction and radia- 
tion). It Is anticipated that any future modeling effort will rely heavily 
on that work. 

While all the existing models so far have assumed thermodynamic equi- 
librium conditions, it might be necessary to consider nonequilibrium condi- 
tions in the futute, as pointed out by Batteh and Keefer (Ref. 4) in their 
discussion of the discrepancy between observed and predicted absorption wave 
propagation velocities. A calculation of radiative properties of nonequi- 
librium hydrogen plasma has been presented by park (Ref. 10) in terms of a 
computer code which calculates emission and absorption coefficients from 
given electron temperature, electron density, neutral particle density and 
intensities of incident radiation. In calculating radiative transport the 
program shows chat there is a large difference between calculated intens- 
ities of radiation emitted by a bulk of equilibrium and nonequilibrium 
hydrogen plasma. 
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3. THRUSTER MODEL REQUIREMENTS 


B«sed on the preceding diacueaion of exiating modela it ia now poaaible 
to eatabliah the baaic featurea which an improved laaer-heated thruater 
model ahould incorporate. Theae featurea are: 

1. Axiaymmetric , two-dimenaional flow ao that effecta of radial 
velocity componenta can he aiaeaaed 

2. Axial aa well aa radial heat tranafer, both by thermal 
conduction and radiation, and 

3. Converging (or focuaed) laaer beam configuration to 
determine the required maaa flow more accurately. 

Requirement (2) haa implicationa which need further elaborationa. It will 
be ahown that thia leada to additional requirementa. 

Aa previoualy diacuaaed, the propagation of plaama radiation (VUV) 
upatream to the front of the abaorption wave will likely lead to a non- 
equilibrium ionization aituation. Itie reaaon ia that electrone have the 
apecial property that their particle maaa ia much leaa that of any of ^he 
other conat ituenta. Becauae of the inefficiency of energy exchange in 
elastic collisions between particles of disparate mass, the electron temper- 
ature may differ appreciably from the temperature of the heavy particle 
species (Ref. 11). Additional equations must therefore be included in the 
model to account for electron mass and energy conservation. 

A second, and probably more serious implication of requirement (2) is 
the fact that the inclusion of upstream influences renders the mathematical 
problem elliptic. This means that a Molution cannot be attempted using a 
single pass forward marching scheme. Instead, the problem haa to be solved 
considering the entire domain in a time-dependent fashion. 
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Ttie ataced requirenenCa render Che potential nodel acre difficult Co 
evaluate chan previou' nodela, particularly ao if it ia viaualized Chat Che 
baaic aet of equationa governing the fluid acchanica (conaervation of maaa 
and energy plua Che Navier-Stokea equationa for conservation of moaentum) 
will have to be aolved numerically. 
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4. MODEL FORMULATION. ANALYSIS AND RESULTS 
4.1 DETAILED MODEL FORMULATION 

A major goal of this effort was to investigate the effect of two- 
dimensional fluid mechanics on the conditions in a laser heated thruster by 
utilising the available l/>ckheed CIM code (Ref. 12) which numerically solves 
the full Navier-Stokes equations (in conjunction with mass and energy con- 
servation equations). Using this code we can therefore satisfy requirement 
1. and. as far as thermal conduction is concerned, also requircownt 2. as 
stated in the previous section. 

The basic equations for axisymmetric flow as used in th« CIM code arc: 


£U ♦ ^ 

5t ax r ar 


(rF) 


♦ C ■ 0 


( 1 ) 


where U, E, F and C each represent four components for the equations of con- 
versation of mass, momentum in x- and r-directions and anergy aa follows: 


U = 


P 

pu 

pv 

p€ 


E 


pu 

2 

pu ♦ p - T 

XX 

puv - T 
^ xr 

u(p^ ♦ p) - uT^^ 


vT 


xr 



( 2 ) 


(3) 
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where 


pv 


F = 


puv . r 


xr 


rr 


pv 


V p) - UT _ - VT _ - k 


xr 


rr 


ST 


G = 


0 

0 

^ 4. ^69 
dr r 


( 4 ) 


(5) 


e - 

e ♦ i (u^+ v^) 

(6) 

'^xx ■- 

2M AP.V 

QY^ 

(7) 

‘»’rr ■- 

dv ^ 

2M ~ ♦ xr * V 

(8) 

^xr ~ 

u ( in ♦ iX) 
^ ' Sf ax' 

(9) 

6* 

0 

II 

2P ( y ) f xr • ^ 

(10) 

V’V z 

3u . 1 d , . 

rtl* 7 ■ 57 ''''' 

(11) 


The above equations represent four equations for the four unknowns - p, u, v 
and T ' if we express the pressure through the equation of state as a func- 
tion of p and T, i >e. , 


P - PRT 0?) 

which is vslid for a perfect gas. The tern Q in the energy equation (see 
Eq.(^)) denotes the net rate of energy addition to the flow which is repre- 
sented by the difference between laser energy absorbed and plasma radiation 
emitted. Since a primary goal was to compare out results with those 
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obtained from Che Batteh and Keefer analysia (Ref. 4), we utfe their 
formulation to eapreae Q aa 

Q - ai^ e"^^- m k (T - T^) (13) 

wliere repreaeiiCa the incident laser energy flux of a collimated beam 
and Q is the absorption coefficient. The second term represents a rather 
crude approximation for the energy loss due to plaama radiation, with the 
form of this term chosen so as to make the Batteh-Keefer anayCical solution 
possible (op. ci*..). 


AC Che present time, the CIM code is formulated in terms of constant 
Chcnsodynamic and transport properties. The rather wide range of tempera- 
tures of interest in the lase.' propulsion problem therefore requires us to 
select suitably averaged quantities for the specific heat, the thermal 
conductivity and the viscosity. In particular, since we want to compare 
results with those of Batteh and Keefer, we want to use the same average 
vjiliit* tor the ratio of sprcitic heat to thermal conductivity that was used 
in their calculations. For consistency, the same averaging procedure is 
applied to Che gas constant. 


Figures 3 through 5 show the thermal conductivity (Ref. 3). the 
specific heat at constant pressure (Ref. 13) and the owlecular weight (Ref. 
13). respectively, for equilibrium hydrogen at p ■ I atm as a function of 
temperature. Integrated averages for c^. k and the gas constant R are 
also shown. From these. w« • in obtain, for use in the CIM code calcula- 
tions, average values for the ratio of specific heats and the viscosity by 
applying standard relationships such as 


and 


y . 


_2 

R 
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Thermal Conductivity. 
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Flit. T - Hy<i roiicn Thermal Cumluclivity (for equiltbruim at 
p 1 atm, Hef. 5 ) 
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with Che latter taken from Svehla (Ref. 14). Finally, Che constant m in the 
radiation loss function has been evaluated by Keefer eC al., to be 

_2 

o ■ 14.3 ca (at p ■ 1 atm). 

Before we atteapc a numerical solution of Che full set of equations 
given in this section, it is useful to investigate possible siapli'.ications. 
As it turns out, for the conteaplated conditions such simpiiticat .ms are 
possible. With additional assumptions, this leads to an e >iicic analytical 
solution which can easily be evaluated. This will be discussed in the 
following section. 

4.2 A SIMPLE ANALYTICAL SOLUTION 

While it is our primary goal to obtain a nuawrical solution of the full 
set of equations presented in the previous section, obtaining such a solu- 
tion will be facilitated by some advance knowledge of what that solution 
might look like. With this in mind, let us non-dimnsionalise Che equations 
by introducing suittible reference quantities for all relevant variables. 
Drfining dimensionlfss variables as 


P’ 

= P/Po 

u', v' 

^ u. v/u^ 

P*. A’ 

- P. A/P^ 


= X, y/t, 1' 

T* 

= t/t^. 

Q * 

V 


k* 

= k/k 


o 


■ V' 



T 

o 


(16) 


and 

for 


introducing these *nto Che equations, Che equations 
convenience) take on the following form; 


(with primes dropped 
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Continuity : 


if ♦ -g^ (P'1) > - 37 J- (npv) = o 


x-Momentum: 


4 <pu, , 4 fpu^ 'j 


r- Momentum: 


X •’ r 

T ^ I 



(IB) 

■ IT. %n] 




1 , 


n 00 ” 

(l'») 


TIT ' 7i7 [‘* ^ Ep) 

I?>7 "If] ♦ |'>[' 

-■fe TilFr"l7]j 


V Ep) 


p ^ 

V 


o 


(20) 


where we have defined 
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Re 




(Reynolds number) 


and 


Pr 


E 


M c 

'^O VO 



Ihie to Che non-dimensxonaliaing proceas all flow variable groups are 
now of order unity, and we muat investigate Che order of nagnitude of the 
dimensionless parameters Re, Pr and E. While Pr is generally of order unity 
for all gases, it can be shown that, for the conditions in the laser heated 
thruster, the Reynolds number is also of order unity, while E is several 
orders of magnitude smaller. We can conclude chat all terms are of equal 
importance in Che continuity and momentum equations, but in Che energy 
equation, however, the pressure terms, all carrying Che E as a multiplier, 
will be orders of magnitude smaller than the convection and conduction terms. 


In pursuit of a simple analytic solution we now examine steady state, 
inviscid one-dimens ions 1 flow at constant pressure. With these assumptions, 
the governing equations reduce to 


Continuity: 



0 


( 21 ) 






Q 


( 22 ) 


Recognizing thee (. t * p/P ) represent* Che total enthalpy which, for our 
case, consists of mostly thermal energy, we ignore the kinetic energy and 


= H=mh=cpT + constant 


Knowing that Che CIH code presently requires c^ and k to be constant, the 
energy equation reduces to 

lii dT _ Q 

■ IT” dx * * T 


where, from Eq. (13) 


Q = ale' 
o 


mk (T.TJ 


Of course, this ordinary differencial equation for T(x) is almost the same 
equati.'tn used by Raiser (op, cit.) except tnat we further simplified it by 
introducing an integrated average k and by neglecting radial conduction. As 
a matter of fact, we can easily include radial conduction in the same manner 
as Raizer did by modifying the constant m in the plasma radiation loss term. 

We can rewrite the differential equation at 


d^T dT _ 

— r - a -r- - m T = 

, Z dx 

dx 


P -ax _ 
e - m T 
k 1 


and arbitrarily specify that the laser power density P ■ obeys 

P =0 for X < 0 


= P for X > 0 
o — 


i- 
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such that X ■ 0 designates the location of plasma ignition (onset of ioniza- 
tion and laser energy absorption). Solutions are then obtained separately 
for X < 0 and x > 0, and the solutions and ?2 joined by requiring 
that both T and its derivative match at x ■* 0. Thus, the boundary condi- 
tions take the following form: 


X = 0: 



x 


JO 


X — ♦ J5 



1 

0 


Standard procedures for nonhomogeneous linesr differential equations then 
yield for x ^ 0: 




and for x > 0: 


af b 


T, (x) = (T - T ) c 


♦ T 


u 1 


(27) 


T^lx) 


1^’ 


\ T. 

d 1 


(28) 


where 



b ♦ a Za 
b - a - Za 

(To-T.) 

(29) 

c 

Zb 




b • a • 2a 

(30) 

b m 

^ + 4 m 

> 0 

(31) 


Imposition of the boundary conditions on T^Cx) generates an additional 
relationship for the laser power density as a function of a, b and a, i*e., 


P 

o 


jk [m -a(g<. a)l b 

b - (a f 2a) ' o‘”i' 


(32) 


The above equations show that the aolution has one free parameter, namely 
the "ignition" temperature T^. Alternatively we could consider the 
initial laser power density as given and determine the ignition 
temperature required to satisfy the boundary conditions. 

Results for typical conditions represented by a, a, and a as used by 
Keefer (Ref. 15) are shown in Fig. 6. Tbe initial temperature, T^, was 
.ifftiinM'tl I <1 !*«■ U)(l K, (iiitl (he ignition tem|>«'toluro was taken ns 9U00 K 
representing the onset of ionization for hydrogen at p •• 1 atm (see Fig. 

5). The teupersture distribution through the laser absorptio- wave is in 
excellent qualitative agreement with results obtained by Keefer (Ref. 15) 
for the ssu.e choice of parameterb. Of course, since our results represent a 
strictly one-dimensional situation (i.e., radial conduction and finite laser 
beam dimensions are ignored), exact quantitative agreement with Keefer's 
results cannot be exoected. 

Another point learned trora this one-dimensional exercise is that, with- 
out the explicit inclusion ui the radiation loss term in the differential 
equation, the differential equation does not have a solution which is able 
to satisfy the downstream boundary conditions. Both the temperature and its 
derivative will become infinite as x approaches infinity, i.e., the soluiion 
diverges. Even though the full conservation equations used in the CIM code 
contain the radial conduction term, and therefore provide an energy loss 
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Fii». 6 - Temperacure Distribution Through Absorption Wave (Hydrogen at p = 1 atm) 


Icrni, (liiti m.iy iiol bo ol Hullicient aiugnilutlu lo prcvciil thu tent|u:rdLiiro 
Iron reaching unreasonable values. As can be seen Iron the example 
presented by Keefer, the power loss due to plasma radiation far outweighed 
that due to thermal conduction in the radial direction. 

The purpose of investigating this simplified analytical model, although 
it duplicates to some degree work previously accomplished, was to provide 
insights which should be of help to obtain a solution from the CIM code. 

This will be discussed in the following section. 

/«.! TWO-DIMtNSIONAI. FLUID MECHANICS EFFECTS 

The successful use of detailed numerical methods such as those used in 
the CIM code requires caretul planning with regard to grid selection, 
initial values and boundary conditions, time step selection and possibly 
scaling considerations to accelerate convergence. 

The analytical results shown in the previous section (as well as 
Keefer's results) show that we have to deal with very steep temperature 
gradients in the laser sbsorp'^ion wave. This imnediately suggests a rather 
tight spatial grid, but how tight? Selecting too coarse a grid can generate 
oscillations in the solution; but too fine a grid will be costly it. terms of 
computation time and storage. 

Using a t imf-dcpcndent metiuid tor the coroputal ion of subsonic Ilow 
fit* Ids, particular c.iro must be exercised in order to ensure that boiindnrv 
conditions remain well iletincd during the course ot the calculation. Ihe 
laser heated thruster is modeled as a straight circular duct with vanishing 
velocity at the wall, spccitied uniform velocity at the inflow boondaiy and 
an outflow velocity distribution to be determined by the analysis. In order 
to keep the problem defined, the laser absorption wave must be positioned 
far enough downstream from the subsonic inflow boundary such that the latter 
IS nut subjected to changes caused by the upstream influence exerted by the 
former. The analytic solution discussed previously shows that, for the 
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sample conditions chosen, the laser absorption wave must be positioned in 
the duct no closer than approximately two cm to the inflow boundary (see 
Fig. 6), assuming that the analytical solution ia quantitatively correct 
with respect to gradients expected in the axial direction. 

Initial conditions, and average values for thermodynamic and transport 
properties must be chosen carefully to be representative of the enormous 
temperature range that has to be covered (with constant values) and, simul- 
taneously, they must be consistent in representing the parameters used in 
the calculations with which we want to compare the CIH code results, vis. 



5.0 cm 


-1 


-f 


- 2.75 ♦ 10^ cm aec/g 


(33) 

(34) 


A consistent set of variables can be selected as follows: According to 
Keefer's results wo expect a maximum temperature of around 20 *10^ K. 

Frt>m Fig. 4 it is found that an integrated average value for the specitic 

heat covering this temperature range ia c~ ■ HO J/gK. Tbe given ratio of 

— — -2 ^ 

c /k then yields k ■ 4.0 • 10 J/cm sec K, a reasonable value as can be 
P 

seen from Fig. 3. While the analytic solutions are determined on the basis 
of global parameters such as pu c^/k, the CIM code requires specification 
of these variables separately. In fact, the Girt code uses y and R to com- 
pute the specific heat. Selecting an integrated average of R ■ 9.4 J/gK 

from Fig. 5, we can compute the ratio of specific heats to be 1.094. The 
equation of state and Eq. (34) then yield 3.593 *10 g/cm and u • 
50.602 cm/sec, assuming p * I atm and T ■ 300 K as initial conditions. 


The selection of the spatial grid for this problem (21 stations in the 
axial (iirection, 26 points in the radial direction) represents a compromise 
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between anticip«te<l requiTement* concerning rctolution end the limitation* 
of Che HSFC univac 1100 *y*te« with respect to compuCacional speed and data 
storage. 

Pure flow calculations without laser energy addition were perfo«^med 

fust. Results from these can be checked against exact (i.e., analytical) 

solutions for axisymmetr ic Poiseuille flow. Initial calculations with a 

time increment as Jictated by the CFL criterion showed extremely slow 

convergence. Huch faster convergence was eventually realised by scaling the 

problem. Scaling as applied here involves multiplication of key variables, 

i.e., the velocity, the viscosity and Che thermal conductivity in this case, 

1 

by a scale factor (1 • 10 , typically) such that Reynolds number and 
Prandtl number are preserved. This type of scaling does increase the Mach 
number. However, because of Che rather small velocity, even the scaled Hach 
number remains saich smaller chan unity, and therefor* th* basic subsonic 
flow character of the problem ia preserved. 

Typical results sre shown in Figs. 7 and d. As seen in Fig. 7, vhc 
initial velocity profile is uniform (over the entire spatial domain) except 
for the velocity at the wall which is set Co sero. As Che parabolic velo- 
city profile developes in time slong the channel, its magnitude in Che 
out flow plane is continually adjusted by mass balance considcratio** '. This 
.id luittmenc causes a feed-back into Che interior of the flow field a.»J 
(hereby enhances convergence of (he solution. Figure 7 shows Che oi^tflow 
velocity profile to have the prop^'r parabolic distribution with a maximum 
velocity on the centerline at roughly twice the value of the average velo- 
city (as represented by the inflow velocity distribution) as predectad by 
the exact solution. A complete velocity vector plot is presented in Fig. 8 
which very clearly shows (he contraction and acceleration of Che tlow in the 
interior of the channel due to the boundary layer uevfloping along the 
channel wall. 

These results also demonstrate well the viscous nature of the flow. 

Note that the scaling factor applied does not affect the balance between 




Fig. 7 - Axi&l Velocity Distribution in Axisymnnetric Flow Charnel (without laser heating) 
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VISCOUS and inurtial furcaa because the Reynolds number is preserved. Wacu 
iiilrutiuciiiK lliv Hcalinx tactur into the equations it can be seen that its 
fllfcl la to increase the viscous dissipation terms relative to the convec- 
tion and conduction terms in the energy equation. This gives the scaling 
factor basically the function of a convergence accelerator, particularly so 
since it was shown previously that the dissipation terms are negligible for 
very slow flows, and remain so in spite of the scale factor. Typical calcu- 
lations involving a time step of lit ■ 1 * 10~^ *cc, and a thousand time 
steps to obtain a converged solution require about one hour of Univsc 1100 
computer time. 

All efforts to obtain a numerical solution for the flow with laser 
heating were unsuccessful. In order to gam an understanding why solutions 
could not be obtained, it is useful to take a closer look at the time step 
criteria which must be observed to control the stability of the numerical 
calculations. The first one is the CFL condition, which can be expressed as 

^*CFL ^ “|u a (35) 


where ^ is the smallest spatisl grid spacing used, and s is the speed of 
sound. The second is a thermal stability criterion, which, when specialised 
to flows with constant pressure (which is very nearly true for our case), 
reads 


.lit 


T 





-I 


h 



(36) 


i 


specialized to our case for constant properties and grid spacings of “ 
0.2 ca and ^r « 0.05 ca, these criteria bccoae 


^‘CFL - 


2.25 • 10 


-7 



(sec) 


(37) 


At 


T 


< 


1.16-10“^ 



(38) 


It is iamediately evident that i!)t^ presents the store stringent criterion, 
especially for flows with high teaperaturcs as they are expected in th« 
laser heated propulsion problem. Since typical cases, using Jt •• 1 * 10~^ 
(sec) require about one hour to converge on the Univac 1100, a teaperature 
rise to about 20,000 K will reduce the peraissible time step by nearly two 
orders of magnitude. Assuming that the total time required to obtain a 
converged solution remains the same as before, we are faced with computer 
times of the order of a hundred hours per case. This is clearly an impos- 
sible proposition. The only remedy here appears to be a computer such as 
Che STAR or Che CRAY, both of which can better handle the storage require- 
ments imposed by this problem, and both of which work at much faster compu- 
tational speed. Additional computational time savings can most likely be 
realized by utilizing a variable time step Co be evaluated as a function of 
temperature as given by the time step criteria. This feature is included in 
the CIM code version for Che STAR computer. Unfortunately, the latter (with 
Its CIM code version tailored to it) was not available for this study. 
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5. CONCLUSIONS 


The numerical analysis of two-dimensional fluid mechanics effects in 
the laser heated thruster is characterized by low density, very high temper- 
atures and extremely large gradients in both variables. Because of subsonic 
(low throughout, the governing equations are of the elliptic type and 
require a solution by a time relaxation method with extremely small step 
sizes in time as dictated by the stability criteria. The first prerequisite 
for a successful solution of this problem therefore is a computer facility 
which can handle the storage requirements and is sufficiently fast so that 
the problem can be solved in reasonable time. The Univac 1100 system does 
not meet these requirements. 

Simplified analytical solutions have shown that there exists a partic- 
ular relationship between the gas flow velocity and the laser power which 
the flow can absorb and simultaneously satisfy the downstream boundary 
condition of a finite temperature. While this relationship is easily 
obtained as part of the analytical solution of the simplified equations, 
most likely it can only be satisfied in a numerical solution via an 
Iterative approach. This implies repetitive calculations and therefore 
amplifies the requirements for a large and fast computer. 

Numerical solutions of subsonic flow problems via time dependent 
methods require particular attention to the specification of boundary 
conditions at the inflow and outflow boundaries. This is a subject of 
current research. It is not known at present to what extent the lack of 
appropriate boundary conditions has contributed to the failure of the 
present effort. 

The fact that the CIH code version usable on the Univac IlOO syst'^is is 
restricted to constant thermodynamic and transport properties should by 
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itsi'lf not prevt'nt a solution. Clearly though, when the temperature vanes 
to the extent as it does in the laser heated thruster problem, the use of 
temperature dependent properties would appear to be more realistic. 

Finally, considering the critical dependence of the permissible time 
step on the temperature, any future attempts at solving the flow problems in 
the laser heated thruster should use a variable time step methodology in 
order to minimize computer time (and cost). 
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